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Abstract 



The theory of the ultrasound propagation in the mixed state of type-II super- 
conductors is suggested which takes into account the Magnus force on vortices, 
the anti-Magnus force on ions, and diamagnetism of the mixed state. The 
acoustic Faraday effect (rotation of polarization of the transverse ultrasonic 
wave propagating along vortices) is shown to be linear in the Magnus force 
in any regime of the flux flow for wavelengths used in the ultrasound experi- 
ments now. Therefore, in contrast to previous predictions, the Faraday effect 
should be looked for only in clean superconductors with a large amplitude of 
the Magnus force. 

Investigation of ultrasound propagation in the mixed state of a type-II superconductor 
has proved to be an effective method of studying high-Tc superconductors (see also Ref. 
for ultrasound experiment and theory in low-Tc superconductors). Recently Dominguez 
et al. have attracted attention to an interesting manifestation of interaction between 
the ultrasound and vortices: acoustic Faraday effect. The effect results from a force on 
ions transverse to the ion velocity. Due to this force the velocity of the transverse sound 
propagation depends on the sign of circular polarization and the plane of linear polarization 
should rotate when the sound wave propagates along vortices. 

However, the system of equations used by Dominguez et al. (see also their later paper 
^) is not Galilean-invariant and does not satisfy the momentum-conservation law for the 
whole system "ions+electrons" since they missed to include into the ion equation of motion 
the so-called anti-Magnus force which is especially important for dirty superconductors. 
In the present work I use dynamic equations which are free from these deficiencies. The 
correct theory predicts other functional dependences on the physical parameters and different 
conditions for observation of the Faraday effect. In particular, our analysis does not confirm 
the conclusion of Dominguez et al. on the strong Faraday effect without the Magnus force 
in dirty superconductors at high temperatures and predicts a many orders of magnitudes 
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weaker effect at low temperatures. Finally, for the ultrasound wavelengths used in the 
experiments the best conditions for observation of the acoustic Faraday effect are at high 
temperatures in clean superconductors. The present work derives the theory of Refs. [jl],|^ 
as a particular limiting case, but generalizes the theory on a much larger domain of physical 
parameters: the nonzero Magnus force (clean superconductors), low magnetic fields, where 
the diamagnetism of the mixed state is essential, wavelengths short compared to the London 
penetration depth. 

We adopt the following picture of electron and ion motion induced by the ultrasound. 
As usual, the electron liquid consists of two parts: normal and superfluid. The normal 
component is effectively clamped to the crystal ions by viscous forces responsible for the 
normal resistance: they move together with the ion velocity Vi = dui/dt {ui is the ion 
displacement). In contrast, the superfluid electrons move with the superfluid velocity Vs 
which is different from the ion velocity in general. Such a physical picture holds until the 
magnetic field is weak compared to the upper critical magnetic field ifc2 and the fiux-fiow 
resistance is much less than the normal resistance. Then the normal current proportional 
to the velocity of the normal electrons with respect to ions may be neglected. 

Thus our three-component system (ions, normal and superfluid electrons) becomes effec- 
tively two-component as in the two-fluid model for superfluids: there is a superfluid with the 
velocity Vs and the mass density meUs and a heavy normal fluid with the velocity Vi and the 
mass density min + me{n — ns) ^ rriin. The charge densities of these two components are eUg 
and —eUs- Here e is the electron charge, rrif. and rrij are the electron and the ion masses, and 
Us and n are the superfluid electron and the total charge number densities respectively. We 
can write the equations of motion for the superfluid and the normal fluid immediately using 
a close analogy with the two-fluid hydrodynamics for rotating superfluids (modifications due 
to presence of the electromagnetic forces are self-evident) 0]: 
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dt'^ * " rrii n V c ' " 

+ -^[{vL-v,i)xB]^ , (2) 

Here Cf is the "bare" sound velocity ignoring interaction with superfluid electrons, the equi- 
librium magnetic induction B is proportional to the vortex density B/^q, $o is the magnetic 
flux quantum, = dui/dt is the vortex velocity, ul is the vortex displacement, and Vsi is 
the local superfluid velocity determined at the points of the vortex line which differs from the 
average superfluid velocity Vg because of deformations of the vortex lattice. If the transverse 
sound wave propagates along vortices (the axis i), then 



Vsl = Vs-\ ^^(^44 



z X 
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The renormalized tilt-modulus CI4 relates only to the vortex line-tension, without in- 
cluding the elastic energy of the average magnetic field ||^ . For an isotropic superconductor 
not very close to Hd Cl^ = (<l'o-B/47rA^) ln(a/rc) (for an anisotropic case see, e.g., and 



references therein). Here a ~ y^o/B and Tc are the intervortex distance and the vortex 
core radius. 

The electrical field E and the magnetic field h generated by the ultrasound wave satisfy 
the Maxwell equations: 

where j = ens{vs — Vi) is the average electric current. The superfluid component of the 
current, the supercurrent js = engVs, should satisfy the London equation averaged over the 
vortex array cell: 

47rA^-± ^ 

V X j, = b^-h . (5) 



Here A = ^ rrieC^ /ATie'^ris is the London penetration depth and is the a.c. component 
of the vortex field (the vortex induction in 0) which coincides with the magnetic field 
h only in an uniform vortex lattice. For the transverse sound propagating along vortices 
6„ = Bdui/dz. 



We need also the equation of vortex motion which connects three velocities Vi, Vsi, and 

Vl- 

- r]{vL - Vi) + r]'[z X {vl - Vi)] 

= 7Thns[z X {Vsl - Vi)] . (6) 

Equation (P) is a self-evident generalization of the equation of the vortex motion known for 
a crystal at rest: the velocities vl and Vgi are replaced by the relative velocities vl — Vi and 
Vsl ~ Vi in accordance with Galilean invariance. The Lorentz force on the right-hand side is 
balanced by the viscous force oc rj and the Magnus force oc rj' . The parameter rj' was known 



to vary from nfins for superclean superconductors |10| to for dirty superconductors [11 
However, rj' may even exceed tt/ots, being equal vr^n, when the lordanskii force is essential 
(see discussion and references in Sec. X of Ref. [^). 

Equations ([l|)-(||) is a closed system of equations averaged over the vortex-array cell. 
Equation (|l]) for electrons does not differ from the case of the ion lattice at rest. As for Eq. 
(0) for ions, I shall show now that the force from vortices in this equation (the last term on 
the right-hand side) is required by the momentum conservation and the third Newton law. 

If there is no external forces on the electron superfluid, the Helmholtz theorem tells 
that Vl = Vsl HI- This means that ?7 = and rj' = tiUhs- Bearing this in mind, one can 
present vl as vl = Vsi + {vl — Vsi) in the right-hand side of Eq. (D). Then the contribution oc 
{{vl — Vsi) X B] is an external force due to interaction with the normal fluid. Correspondingly, 
the same force must appear in the ion (normal fluid) equation (|]), but with the opposite 
sign. In superfluid hydrodynamics this force is called mutual friction force But for 
superconductors it is better to call it coupling force, since the force may incorporate not 
only friction, but also elastic pinning (see below). The coupling force, being external for the 
superfluid component, is internal for the system "ions-|-electrons" as a whole. Therefore this 
force doesnot contribute to the time variation of the total momentum of ions and electrons 

— [miUVi + me[n - ns)Vi + men^f^J 



Here ji = ens{vsi — Vi) is the local current on the vortex hne, in contrast to the average 
current j = eusivg — Vi) in the Maxwell equations (^). Equation demonstrates the 
momentum conservation law for the case when the wave propagates along vortices (the axis 
z): the time variation of the momentum is determined by the divergence of the momentum 
flux which consists of (i) the stress tensor of the ion lattice, (ii) the stress tensor of the 
magnetic field linearized with respect to the a.c. field h, and (iii) the vortex-lattice stress 
tensor given by the tilt-modulus C44. In the long- wavelength limit AA; ^ (see below) the 
superfluid component oc dvg/dt of the momentum variation may be neglected, and Eq. (^ 
has a simpler form: 



{d''ui/dt - clV^u^ = X 5] . (8) 



This equation was used in the analysis of Shapira and Neuringer |0, but they neglected 
difference between the local and average electric current which is responsible for the dia- 
magnetism of the mixed state. 

Replacing the coupling force in the equation of ion motion Eq. (^ with the help of Eq. 
(I) one obtains: 

d'^Uj r,-Lr,^ e n, f -± 1 , 



nii n [ c 



-^VivL - Vi) - ^^^^[{VL -Vi)xB]^ , (9) 

where fj = rj/nhus and om = ri'/Trhus are the dimensionless amplitudes of the viscous and 
the Magnus force respectively. The coupling force has a component oc [om — 1) transverse to 
the ion velocity which is maximal in a dirty superconductor when the Magnus force vanishes 
{ttM = 0). Indeed, suppression of the Magnus force in the dirty superconductor means that 
impurities produce a force which cancels the Magnus force on vortices, i.e. on superfluid 
electrons. Then in accordance with the third Newton law one must expect that the force of 
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opposite sign acts from vortices on impurities, i.e., on crystal ions. This explains the force 
oc — 1) which may be called the anti-Magnus force. 

Later on we assume that the transverse ultrasound plane wave oc exp{ikz — iujt) propa- 
gates along the vortices. From Eqs. (0)-(|^) one can derive relations connecting the electric 
field and the local superfluid velocity with the ion and the vortex velocities and displace- 
ments: 



m, 



^ _ _ ■■^e _2 



(10) 



Vsl 



1 + erisB 



Cu{k)k'^[z X ul] . 



(11) 



where now 



C44(fc) 



S2 



An 1 + A2A;2 



+ c: 



44 



(12) 



is the fc-dependent tilt-modulus related to the total energy of deformation |T^. The contri- 
bution oc uj"^ to the electric field E will be neglected later on, since it is not connected with 
vortices. It yields a small correction to the sound velocity of the relative order niens/min 
which is present even in the Meissner state. 

Finally, with the help of Eqs. (^, (|l3), and ([TT|), one can rewrite the equation of the ion 
motion, Eq. (^, in the terms of the ion velocity and displacement only: 

en.B 



[f\\Vi + flXz X Vi]] 



(13) 



The longitudinal and the transverse forces on ions from vortices are given by the dimension- 
less force parameters: 



/ll 



Dk'' 



1 



iuj [t] — Dk"^ / iujY + a\j 



V 



1 + A2A;2 



Dk^ 



21,2 



X'k 



(14) 



2A2p 



+ 7^ 



1 



1 + A2p {j]-Dk'^l%ujf^a 



M 



if De\ x^e Dk 

X 2 77 



+aM 



(15) 



Here D = cC^^/engB. 

The previous experimental and theoretical investigations dealt with the long- 

wavelength case Xk — > 0. One may call it the electrodynamic limit since in this case all 
forces from vortices on ions can be expressed in terms of the electrodynamic parameters: 
the Ohmic and the Hall conductivities, ctq = rjc^/^QB, an = rj'c^/^oB, related to the 
viscous and the Magnus force respectively, and the magnetic permeability 

B^Att _ B^ 
^ ~ By 471 + Q4 ~ 47rC44(0) ' ^^^^ 

which describes the diamagnetism due to circular currents over the vortex-lattice cell 0. 
Here 6*44(0) is the tilt-modulus in the limit A; — (the Labusch modulus). One can check 
that Eq. (|T6|) yields the differential magnetic permeability fi = B/H for the case when the 
variation of the magnetic field is normal to the vortices. Here H = {1/ ATi)dF{B) / dB is the 
thermodynamic magnetic field along the equilibrium magnetization curve. In the terms of 
electrodynamics the equation of vortex motion (P) is simply the Ohm law 

Jz = aoEi -anizx Ei] , (17) 

which connects the local current ji = j/yU and the electric field Ei = E + {1 / c)[vi B] in the 
coordinate frame moving with the ion velocity. Then Eq. (|13D may be derived from Eq. (^ 
together with the Ohm law Eq. (p!7|) and the relation E = {At^iuj / c^k"^)] = {Atiiuj^/ (?k'^)ji 
which follows from the Maxwell equations (^). This yields the following values of the force 
parameters: 



B 



2 



{po — c'^k'^/4:7riujjj,) + (t|^_ 



c^k"^ 

(18) 
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/, = ^f^V . (19) 

One can easily check that Eqs. (|1^) and (|19|) coincide with Eqs. (|T^ and (|1^) if Xk 0. In 
the electrodynamic hmit our theory is vahd in a wider interval of the magnetic fields than in 
the general case of arbitrary Xk. Since we use the equations averaged over the vortex-array 



cell, they hold until the wavelength 27r/k exceeds the intervortex distance a ~ \J^q/B. But 
because of A <C this condition may be satisfied even if a ^ A, i.e., rather close to the 
lower critical field Hd where the diamagnetism is important, i.e., /i essentially less then 
unity. On the other hand, in the electrodynamic approach our assumption in the beginning 
of the paper that the current of normal electrons is negligible is not necessary: one can use 
the conductivity tensor taking into account this current. Then the theory is valid even close 

to Hc2. 

For the sake of simplicity we did not include the elastic pinning force into our analysis 
explicitly, but it is easy to do simply by replacing the viscous coefficient rj hj rj — ap/iu in 
all equations. Here ap is the bulk pinning coefficient which may, in principle, depend on the 
frequency as assumed by Dominguez et al. 0]. Thus our analysis holds for any regime of 
vortex motion, either the non-activated flux flow, or the thermally assisted flux flow (TAFF) 
with flux jumps over pinning barriers. But different regimes of the flux flow correspond to 
different expressions for the conductivities. 

One can obtain the results of Refs. from Eqs. (|18D and ([T9|) neglecting the Magnus 
force and the flnite diamagnetism of the mixed state {an ~ «a/ = 0, /i = 1). Note that 
the magnetic permeability /i of Shapira and Neuringer 0] relates to the atomic magnetism, 
whereas the latter was ignored in the present work and our yU is due to the circular currents 
in the vortex-lattice cell and the difference between the local and average currents as a result 
of them. However, flnally the role of yU in the equations is similar in both cases, as one might 
expect for an electrodynamic theory. 

Our analysis seriously differs from that of Dominguez et al. [^: (i) They missed to take 
into account the anti-Magnus force in the equation of ion motion. This caused violation of 
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the momentum conservation law and the wrong prediction for the Faraday effect in the dirty 
superconductors in which the anti-Magnus force is especially important, (ii) The equation of 
vortex motion used in Ref. contained the laboratory vortex velocity and thereby was not 
Galilean-invariant. This means an assumption that the laboratory frame is preferential. But 
for our problem the only preferential coordinate system is that moves with the ion velocity. 

Experiment and theory [||] have shown that in a Bi superconductor there is a crossover be- 
tween two temperature regions: the low-temperature region of high Ohmic conductivity due 
to high activation barriers in the TAFF model where ufj/Dk'^ = Ctfj/Dk = ATcufiao/c'^k'^ > 
1, and high-temperatures region of low conductivity where Anufiao/c^k^ < 1- The predic- 
tions of Dominguez et al. for the Faraday effect must be revised both in the low-temerature 
and the high-temperature regions. At low temperatures they obtained that the Faraday 
rotation (the angle of polarization rotation per unit length) is dO/dz = aHB^/2minctc'^, 
whereas Eq. (|19D yields the Faraday rotation dO/dz = ensBR.ef±/2minctC which is by the 
factor [c^k"^ / Atiuj ^aoY smaller. This factor is of order unity at the crossover temperature 
T ~ 60 K, but at low temperatures the Ohmic conductivity increases proportionally to 
exp([//T) where the activation barrier is about 500 K for typical magnetic fields |T|. Thus 
Dominguez et al. overestimated the Faraday effect at low temperatures by many orders 
of magnitude. At higher temperatures Dominguez et al. obtained the Faraday rotation 
d9/dz = eB/2miCtC which did not depend on the Hall conductivity. This means that the 
Faraday effect is possible in a dirty superconductor without the Magnus force even in the 
limit A/c ^ in disagreement with Refs. Dominguez et al. explained it by the effect 

of the electromagnetic force oc (^E + M^vi x 5]^ [see Eq. (|^) for E] in Eq. @, which was 
neglected in Ref. [0]. However, as explained after Eq. (|), this equation should include also 
the anti-Magnus force which cancels the electrodynamic force in the limit A/c ^ restoring 
the result of Refs. According to Eq. (|T5]) the Faraday rotation without the Magnus force 
(om ~ ct// = 0) is given by 



dO eUsB 



Reh 



eUsB A^P 1 



(20) 



dz 2minctC 



rriinctc 1 + X^k"^ 1 + {c^k"^ j AnojuaoY 
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For the ultrasound frequencies 3 MHz from Ref. [0 this expression yields at high temper- 
atures > 60 K the Faraday rotation at least by the factor 10~^ smaller than predicted 
by Dominguez et al. 0. So the Faraday effect in dirty superconductors is possible only at 
AA; > 1. However, the frequency of ultrasound must be very high in order to achieve this case 
(about a few GHz). For the frequencies used in the ultrasound experiments now, the best 
conditions for observation of the Faraday effect are clean superconductors with large Hall 
angle at high temperatures. Then the Faraday rotation can be strongly amplified close to 
the resonances with vortex modes. The resonances correspond to zeroes of the denominator 



in Eqs. ([T^) and ([191) . Weakly damped vortex modes exist in superclean high-Tc super- 
conductors according to recent experimental and theoretical investigations [[7|,p!3|. Thus the 
ultrasound experiment is able to reveal these vortex modes as has been already discussed in 



Ref. \U 



In summary, the theory of interaction between the ultrasound and the vortices in the 
type-II superconductors has been suggested. The first time the theory takes into account 
the Magnus force on vortices, the anti-Magnus force on ions, and the diamagnetism of the 
mixed state self-consistently. This results in a serious revision of previous predictions con- 
cerning the acoustic Faraday effect. Despite this revision, our analysis confirms that possible 
observation of the acoustic Faraday effect is expected to provide a valuable information on 
vortex dynamics in type-II superconductors. 
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